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The master equation approach to Lyapunov spectra for many-particle systems is applied to non- 
equilibrium thermostatted systems to discuss the conjugate pairing rule. We consider iso-kinetic 
thermostatted systems with a shear flow sustained by an external restriction, in which particle 
interactions are expressed as a Gaussian white randomness. Positive Lyapunov exponents are calcu- 
lated by using the Fokker-Planck equation to describe the tangent vector dynamics. We introduce 
another Fokker-Planck equation to describe the time-reversed tangent vector dynamics, which al- 
lows us to calculate the negative Lyapunov exponents. Using the Lyapunov exponents provided by 
these two Fokker-Planck equations we show the conjugate pairing rule is satisfied for thermostatted 
systems with a shear flow in the thermodynamic limit. We also give an explicit form to connect the 
Lyapunov exponents with the time-correlation of the interaction matrix in a thermostatted system 
with a color field. 
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I. INTRODUCTION 

The Lyapunov exponent is an essential concept to ex- 
press the instability of orbits in a dynamical system. It is 
introduced as an exponential expansion (or contraction) 
rate of an infinitesimal perturbation of orbits, and its 
positivity implies that the system is chaotic. In general 
there is a Lyapunov exponent for each independent di- 
rection of the infinitesimal perturbation of the orbit, and 
the sorted set of such Lyapunov exponents is called the 
Lyapunov spectrum, and has been the subject of study 
in many-particle systems. For example, the existence of 
its thermodynamic limit Jl|-|3| , an effect of the rotational 
degrees of freedom of molecules 0] , its stepwise structure 
and the Lyapunov modes |5|, and a tracer particle effect 
H have been observed and discussed in the Lyapunov 
spectra of many-particle chaotic systems. 

A characteristic of the Lyapunov spectrum that is 
known in Hamiltonian systems, is that the Lyapunov 
exponents appear as a pair, namely any positive Lya- 
punov exponent accompanies a negative Lyapunov expo- 
nent with its opposite sign H. This characteristic, which 
is based on the symplectic structure of the Hamiltonian 
mechanics, is not correct in non-Hamiltonian systems, 
but it is interesting to know how it is modified in quasi- 
Hamiltonian systems such as a Hamiltonian system cou- 
pled to a thermodynamic reservoir. This problem has 
been considered in some thermostatted dynamics with a 
term to extract the heat produced in the system by exter- 
nal force fields, and led to the proposal of the conjugate 
pairing rule for thermostatted systems, which claims that 
the sum of any Lyapunov exponent pair is not zero but a 
constant regardless of the exponent numbers J|,|9| . This 
conjecture was confirmed by many following numerical 
calculations |l(] [l^]. This also led to the discovery that 
some thermostatted systems contain hidden Hamiltonian 
structure |l3|,|l4| • The pairing rule is not only interesting 
as a mathematical structure of the thermostatted system 



but also valuable for a practical use; the conjugate pairing 
rule for the thermostatted system allows us to calculate 
non-equilibrium transport coefficients (eg. conductivity 
and viscosity) from only one pair of the Lyapunov expo- 
nents, such as the larg est and smallest Lyapunov expo- 
nents only ||||Mf. 

A problem is that the necessary and the sufficient con- 
ditions for the conjugate pairing rule to hold for ther- 
mostatted systems is not clearly known. The conju- 
gate pairing rule for the iso-kinetic thermostatted system 
with a color field was proved for the soft core interac- 
tion potential [ and the hard core interaction potential 
fll4] , [l8|] , regardless of the number of particles. A similar 
discussion was done in Nose Hamiltonian dynamics [jl9| . 
These works give the sufficient conditions for the conju- 
gate pairing rule. On the other hand, it was suggested 
numerically that it can be violated in the presence of a 
magnetic field J2£] and in inhomogenously thermostat- 
ted systems such as a system under temperature gradi- 
ent pjj or a system in which the peculiar momenta are 
thermostatted j22|. Another numerical work also sug- 
gested that it is not exact in the iso-energetic thermo- 
stat with a finite number of particles ^3|, although the 
iso-energetic thermostat should be equivalent to the iso- 
kinetic thermostat as the number of the particl es goe s 
to infinity, namely in the thermodynamic limit p4| . p5[ . 
A special interest is the iso-kinetic thermostatted system 
with a shear field, which is described by the Sllod equa- 
tion for the planar Coutte flow |2(J . Early investigations 
supported the c onjuga te pairing rule for such a system 

§, L0|j27| . Refs. p2| , p8| suggested a small deviation from 
e conjugate pairing rule. An analytical consideration 
showed that the deviation from the conjugate pairing rule 
should be at most fourth order in the shear rate in the 
case of a small shear rate in the thermodynamic limit 
p9fl . However a recent numerical calculation with a more 
careful numbering of the Lyapunov exponents and with 
numerical error bars showed that within the numerical 
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precision the conjugate pairing rule was satisfied [ p0[ . Af- 
ter all these trials, a justification of the conjugate pairing 
rule for the iso-kinetic thermostatted system with a shear 
field still remains as an open problem. 

Analytical calculation of the full Lyapunov spectra for 
many-particle systems is still not an easy matter, and 
so far full Lyapunov spectra have been calculated mainly 
using numerical approaches. On the other hand, recently 
an analytical approach to the full Lyapunov spectra was 
proposed for many-particle systems with a random in- 
teraction ]3l| . This approach describes the tangent vec- 
tor dynamics by a master equation, and allows the cal- 
culation of all individual positive Lyapunov exponents 
through the average of the magnitude of the tangent 
vector. It was used to explain the stepwise structure 
of the Lyapunov spectrum for many-particle Hamilto- 
nian systems. However a generalization of this approach 
to non-equilibrium systems, especially to thermostatted 
systems, is not known yet. 

The purpose of this paper is to generalize the master 
equation approach to the Lyapunov spectrum to non- 
equilibrium thermostatted systems, and to discuss the 
conjugate pairing rule of the Lyapunov spectrum for the 
iso-kinetic thermostatted system with a shear field. Es- 
pecially we propose a method to calculate negative Lya- 
punov exponents using a time-reversed master equation. 
In this paper we concentrate on the case where the parti- 
cles interact with a Gaussian white randomness. In this 
case the master equation that describes the tangent vec- 
tor dynamics is attributed to the Fokker-Planck equation. 
We also restrict our consideration in the thermodynamic 
limit, in which the fluctuations of the friction coefficient 
can be neglected [p7| . Under these conditions we show 
that the conjugate pairing rule for the thermostatted sys- 
tem with a shear field given by the Shod equation is sat- 
isfied. As a special case we also discuss briefly an explicit 
form to connect the Lyapunov exponents with the time- 
correlation of the interaction matrix in a thermostatted 
system without a shear field. 



II. ISOKINETIC THERMOSTATTED SYSTEM 
WITH A SHEAR FIELD AND ITS TANGENT 
VECTOR DYNAMICS 



We consider non-equilibrium systems with an iso- 
kinetic thermostat. Our consideration includes the case 
where a shear flow is sustained by an external restric- 
tion, and for simplicity we consider a two-dimensional 
system consisting of N particles with the same mass m. 
We introduce q^(i) = {qi j) (t) , q { y j) (t)) T and p^(t) = 

(Px\t) , py\t)) T as the spatial coordinate vector and the 
momentum vector of the j-particle, respectively, at time t 
with the transpose operation T. (Note that all vectors in 
this paper are introduced as column vectors.) Equations 
for q( J ') (t) and p^ (t) are expressed as [p6[ 



dp^(t) 
dt 



dU(t) 
'0qW)(t) 



[ 7 S 2 + a(t)/ 2 ]p»(t) (2) 



where U(t) is the potential energy as a function of 
3 = 1j 2, • • • , N and t only, and we introduce !B 2 fc 
is the 2k x 2k matrix defined by 
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(3) 



with the k x k identical matrix 1^ and the k x k null 
matrix Ofe. Here 7 is the shear rate as an external pa- 
rameter, namely a constant gradient of the x component 
of the local velocity in the y direction, and a(t) is defined 
by 



a(t) 



Ef =1 P W W T ( 



dU(t) 
9q«> (t) 



72 2 p (i) (*) 



Ej=ilp tf) (*)l s 



(4) 



so that the total kinetic energy is constant in time: 

4££=ilP W (*)l7(2m)]M = 0. Eqs. (0) and (|) are 
called the Shod equation for the planar Coutte flow with 
the isokinetic thermostat, and gives the model for the 
system driven by external fields and/or a shear rate with 
an attached heat reservoir which removes the energy gen- 
erated inside the system and maintains the temperature 
of the system constantly in time. As an example de- 
scribed by Eqs. (Q) and (Q), other than the system with 
a shear field, we may mention the color field system in 
which the system consists of many particles with charges 
of different signs and is driven by an external electric field 

us 

In general, the quantity a(t), which is interpreted as 
the friction coefficient, depends on the coordinates and 
the momenta of the particles, so is variable in time. How- 
ever, it is known that the fluctuation of the quantity a(t) 
is small in a system consisting of many particles ]27j] . (For 
a justification of this point by the kinetic approach see 
Ref. which shows that the quantity a(t) fluctuates 
with the order of l/y/N around a fixed value.) Based on 
this fact, in this paper we consider only the system which 
consists of enough particles so that the friction coefficient 
a(t) in Eq. (||) can be replaced by a fixed constant a. 

For a convenience we represent the 4JV- 
dimensional phase space vector T(t) as a vector 

(q£\t),q£\t), • • • , qi N) (t), q^(t), qf\t), • • • , qf\t), 

p£\t),pl?\t), ■ • • J x N \t), pi 1] (t),pf(t), ■ ■ ■ ,pl N) (W. 
Using this notation and the assumption explained in the 
previous paragraph we obtain the equation 



d5T(t) 
dt 



£{t)5T(t) 



(5) 



for the tangent vector 5T(t). Here the matrix C(t) is 
given by 



dq_^(t) __ 1 
Jt " 



p«(i)+72 2 q UJ W 



C„n(j)| 



(1) 



r ,n= I * hN/m 
L[t) ~ I R(t) * 



(0) 
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with 2N x 2N matrices $, \& and R denned by 

$ = 7 S 2Ar , (7) 



- OtI 2 N, 



(8) 



#(t) = - 



6> 2 E/(t) 

0q(t)0q(t) 



(9) 



where we introduced q(£) as a vector (q£\t), q^\t), ■ 



i N \t),qV(t),qW(t),...,qi N Ht)y 



q.r. 



III. RANDOM INTERACTIONS AND MASTER 
EQUATIONS FOR THE TANGENT VECTOR 
DYNAMICS 

In this section we introduce a random interaction be- 
tween the particles, and obtain the two kinds of mas- 
ter equations corresponding to the time-forward tangent 
vector dynamics and the time-reversed tangent vector dy- 
namics by using Kramers-Moyal expansion technique. 



2N 2N 

EE 

fl—l v—\ 

2N 2N 2N 2N 



dSp^ 

+ E E E E Tpn/twMiAiv 



£1=1 V=\ fl' = l v' = l 



d 2 



dSpfj,>dSp v > 



P M(sr,t) (12) 



applying the Kramers-Moyal expansion technique to the 
dynamics (^|). Here Sqj and Spj are the j-th components 
of the coordinate part 5q and the momentum part Sp in 
the tangent vector ST = (5q, Sp) T , respectively, and <I> M „ 
and \& ^„ are the matrix elements of the matrix <E> and ^ 
defined by Eqs. (0) and (0), respectively. The derivation 
of Eq. (|l2|) is given in Appendix [X|. Eq. (|L^) in the 
special case of $ = 2 n and = 2 n have already been 
used to discuss the stepwise structure of the Lyapunov 
spectrum for a many-particle Hamiltonian system . 



B. Anti-Fokker-Planck equation for the 
time-reversed dynamics of the tangent vector 



A. Fokker-Planck equation for the forward dynamics 
of the tangent vector 

We consider the case that each particle interacts with 
the other particles randomly enough so that the matrix 
R(t) = (Rjk(i)) can be regarded as a Gaussian white 
random matrix satisfying the conditions 



(-Rnifi(*l)-^2f 2 (*2) ' 



^/i2n-lf2«-l(*27l-l)) 



0, (10) 



^ ^ ^Mil V 3\ V 31 ^^33 V 33 ^34 V 3A 

Pi 

■■D b 



'Mj 3 



xS(t h 



— 1 V 3ln — \ V'iln l/ 32n 

tj 2 )S(tj 3 ~ tj 4 ) ■ 



■5{t 3 



tj 2 J (11) 



for any integer n and a 4-th rank constant tensor Djkin, 
where we take the sum over only the permutation Pd : 
(1,2, ••• ,2n) -> (ji, j 2 , ■ ■ ■ , j2n), and the bracket (• • •) 
means the ensemble average over random processes. 

Under the randomness conditions (|l^) and ([□]) the 
time-evolutional equation (J5J) is regarded as a stochas- 
tic equation of the Langevin type, and its correspond- 
ing master equation for the the probability density 
p( + \5T,t) for the tangent vector ST at time t is given 

by 



£,«(«r.O 

2N 2N 



EE 



^a V 8q u 



5p v )p^{8T,t) 



As shown in Ref. pi| , in the case of 7 = and a — 
the Fokker-Planck equation ( |Lq ) provides the positive 
Lyapunov exponents as the time averaged exponential 
rate of the randomness average by the probability density 
p( + ^ (r, t) in the time evolution of infinitesimal perturba- 
tions of the dynamical variables. However this method 
does not provide directly the negative Lyapunov expo- 
nents, because in the stochastic system the randomness 
average of the distance between the infinitesimal nearby 
trajectories should not shrink in the infinite time limit. 
This was not a problem in the Hamiltonian system dis- 
cussed in Ref. |3l|], because in the Hamiltonian system 
the absolute values of the negative Lyapunov exponents 
are the same with the positive Lyapunov exponents M. 
However in the thermostatted system discussed in this 
paper such a simple relation of the negative and the pos- 
itive Lyapunov exponents can not be expected any more. 
In order to overcome this problem and to provide the 
negative Lyapunov exponents using the master equation 
approach directly, we use the fact that the negative Lya- 
punov exponents can be regarded as the positive Lyapunov 
exponents for the time-reversed motion. This fact has 
been used in some works to calculate the negative Lya- 
punov exponents for chaotic systems | fi6| , p4|j35[ . 

In the iso-kinetic thermostatted system with a shear 
field the time-reversed motion is expressed by the " time- 
reversed mapping" t —> — t,q^q, p ^ — p and 
7 — > —7. The transformation 7 — > —7 is justified by the 
fact that the direction of the shear flow changes to the 
opposite direction in the time-reversed motion. This jus- 
tifies the time-reversal transformation a — » —a for the 
friction coefficient by Eq. (^). The time-reversed map- 
ping leads to the time-reversed transformations 
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t - 


-> -i 


(13) 


5q- 




(14) 




-dp 


(15) 


$ - 


> -$ 


(16) 


* - 




(17) 



for the tangent vector dynamics, noting the relations (0) 
and (||) to connect the matrices <1> and ^ with the shear 
rate 7 and the friction coefficient a. It is important to 
note that the tensor Djkin and the tangent dynamical 
equation (y) themselves are invariant under the transfor- 
mations 

Now we introduce the Fokker-Planck equation for the 
probability density p^~ '(ST, t) for the time-reversed tan- 
gent vector at time t as the transformed equation of the 
Fokker-Planck equation ( |l2] ) by the transformations ( p"3| - 
O), namely 



anti- Fokker-Planck equation, but otherwise it can not be 
used to calculate the negative Lyapunov exponents. 

The anti-Fokker-Planck equation is analogous to the 
" anti-Lorentz-Boltzmann equation" which was intro- 
duced to calculate the negative Lyapunov exponents us- 
ing the kinetic approach p^j35| , |36| . In this approach the 
anti-Lorentz-Boltzmann equation is given as the Lorentz- 
Boltzmann equation where the collision operator has the 
opposite sign to the ordinary Lorentz-Boltzmann equa- 
tion in an equilibrium or a non-equilibrium stationary 
state. 

At least, in the equilibrium case expressed by 7 = 
and a — the anti-Fokker-Planck equation must pro- 
vide the negative Lyapunov exponents as the opposites 
of the positive Lyapunov exponents calculated by using 
the Fokker-Planck equation ([l2|). In the next section we 
show it actually as the special case of more general re- 
sults. 



2N 2N „ , _ 



fl— 1 V— 1 

27V 2N 







fl—l v — \ ^ 

2N 2N 2N 2N ^ 



/I— 1 U—l fl' — l u' — l 



d 2 



P (-\sr,t). (is) 



This is simply the equation with the opposite sign of 
the diffusion term to the forward Fokker-Planck equa- 
tion (|l^), and is interpreted as the master equation to 
describe the time-evolution of the tangent vector whose 
initial condition is the time-reversed initial condition to 
the Fokker-Planck equation (fT^). We call this equa- 
tion for the probability density p(~\6T,t) the "anti- 
Fokker-Planck equation" in this paper, and calculate the 
negative Lyapunov exponents as the exponential rate 
of the randomness average by the probability density 
p(~>(5T, t) in the time evolution of infinitesimal pertur- 
bations of the dynamical variables in the minus infinite 
time. 

It is important to note the difference between the anti- 
Fokker-Planck equation ( p"S) a nd the so called "backward 
Fokker-Planck equation'^ 37J . The backward Fokker- 
Planck equation describes the dynamics before an ini- 
tial time in which the initial condition is the same as 
in the forward Fokker-Planck equation to describe the 
dynamics after the initial time. However, in order to cal- 
culate the negative Lyapunov exponents from the time- 
reversed motion, we must use the different initial condi- 
tion which has the opposite sign of the momentum to the 
initial condition in the forward Fokker-Planck equation. 
In the equilibrium case expressed by 7 = and a — 
the backward Fokker-Planck equation coincides with the 



IV. CONJUGATE PAIRING RULE FOR 
THERMOSTATTED SYSTEMS WITH A SHEAR 

We have to know the time evolution of the amplitude 
of the tangent vector in order to calculate the Lyapunov 
exponents. Such a time evolution for the forward move- 
ment (the time-reversed movement) is expressed as the 
dynamics of the diagonal elements of the matrix T (+)(£) 
(the matrix >{t)) given by 



T\(±) 



T (±, (t)H(iq{q T )| 



(19) 



Here the bracket (■ ■ ■) < j ± " > means to take the average by 
the probability density p^(5T, t), namely 



(X(ST)) 



(±) 



d6TX(5T)p^\ST,t) 



(20) 



for any function X(8T) of ST. In Ref. |u| the posi- 
tive Lyapunov exponents were calculated by the time- 
averaged exponential rate of the diagonal elements of an 
orthogonal-transformed matrix of T^(t). We get the 
negative Lyapunov exponents from the matrix T^~'(t) 
by a similar procedure. The introduction of Lyapunov 
exponents by the spatial coordinate part only (or the mo- 
mentum part only) of the tangent vector has been used 
previously by Refs. p5| , |3"8| ] . 

We introduce the matrix T^'^t) defined by 



fW(i) ee T (± )(±i)e 



±at 



(21) 



As shown in Appendix |B], the matrix T^'(t) satisfies the 
differential equation 



d 4 f(±)(t) 



di 4 
1 

+ 16 



dt 2 



dt 2 



n4fO t >(i)-2fi 2 ft | 0(i) (ft 2 



4 



+T^>(t) (ft 4 ) 



2 ^ df(±)(t) 
m 2 I dt 



= (22) 

where we assumed the probability density p^' (ST, t) to 
be zero at the boundary of the tangent space. Here f2 is 
the (27V) x (27V) matrix defined by 

O = $ - (23) 

and T> is the linear operator to map any (27V) x (27V) 
matrix X = (X jk ) to the (27V) x (27V) matrix V{X} = 
( (D{X}) jk ) defined by 



N N 



{V{X}) jk = Y,Y, D i^ X , 



ft is • 



(24) 



It may be noted that Eq. ( p2| ) for the matrix Y^^i) is 
invariant under the transformation Q — > — ft. 

We can choose the initial probability density 
p( + '(<5r,0) arbitrarily to calculate the positive Lya- 
punov exponents. On the other hand in order to derive 
the corresponding negative Lyapunov exponents we as- 
sume the initial probability density p^~\5T, 0) for the 
time-reversed tangent vector to satisfy the condition 
d k f^(t)/dt k \ t =o = d k f ( -+^t)/dt k \ t=Q , k = 0,1,2,3 at 
the initial time t = 0. Under this assumption it follows 
from Eq. (||) that 



T { -\-t)e- &t = T (+) (i)e at , 



(25) 
by 



because the quantities Y' + )(t) and T^'(t) defined 
Eq. ( pl| ) satisfy the same differential equation (22) and 
have the same initial condition. Therefore the diagonal 
element fg°(t) of any orthogonal-transformed matrix of 
T^\t) must satisfy the relation 



T 



(-) 

a 



H) = fW(t)e 2 



(26) 



This equation connects the time-forward evolution and 
the time-reversed evolution in the amplitudes of the spa- 
tial part of the tangent vector in the thermostatted sys- 
tem. 



The j-th positive (or zero) Lyapunov exponent A 

and its conjugate negative (or zero) exponent A 
given by 



(-) 



(+) 
j 

are 



v(±) 



1 f (±; 
lim — In , ]\ , 

i=-»±oo 2t -f (±) 
33 



(t) 



(0) 



(27) 



The quantity Y& J (t) satisfies the condition Y^^O) 



(+)/ 



'33 



. rj (0) at the initial time, and using Eqs. 
we obtain 



33 

and 



4 +, + V 



) 



-a. 



(28) 



This is the conjugate pairing rule of the Lyapunov spec- 
trum for the iso-kinetic thermostatted system with a 
shear field. It is clear that it is attributed to the pairing 
rule of the Hamiltonian system in the case of a — 0. 



V. CONJUGATE PAIRING RULE FOR A COLOR 
FIELD 



For an actual calculation of the Lyapunov spectrum for 
an iso-kinetic thermostatted system by the master equa- 
tion approach we have to know the value of tensor Djkin 
and to solve the differential equation (^2|) for the matrix 
T^\t). Let's discuss these points briefly by using a case 
without a shear field such as a color field system Jl2|,^2| , 
namely 



7 = 0. 



(29) 



For simplicity in this section we also assume that the ten- 
sor -Djfe/n is expressed as the multiplication of the matrix 
elements of a symmetric (27V) x (27V) matrix W = (Wj k ): 



Djkln = W jk Wi n . 



(30) 



This assumption has already been used in Ref. |u| . 

As shown in Appendix under the assumptions ( |29| ) 
and (|30| ) the equation of the matrix TW(t) is simplified 
to 



dt 3 



dt 



WT [± \t)W = {). (31) 



It is shown that Eq. ( |22] ) is given by taking the time- 
differential in both the sides of Eq. (31). Using the or- 
thogonal matrix V diagonalizing the matrix W, namely 



[V T WV) ]k = ujjSjk 



(32) 



with a real eigenvalue u>j, the quantity Yj-^ (t) is ex- 
pressed as the diagonal element of the matrix T^(t) = 
(f^ } (t)) defined by 



Y (±) (i) = V T T {±) (t)V. 



(33) 



Using Eq. ( |33| ) we can solve the equation for the quantity 
fj^\t) derived from Eq. @, and by using Eqs. @, 
( p7| ) and (|33l ) the Lyapunov exponents are given by 



v(±) 



a 1 
— ± - 



A; 



a 
3A 



(34) 



where A,- is defined by 




1/3 



(35) 



(See Appendix |^ about a derivation of Eq. (|34"|).) It is 
clear that the Lyapunov exponents given by Eq. ( |3^ ) 
satisfy the conjugate pairing rule (Eq). 
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Concerning the expression (B4T) for the Lyapunov ex- 
ponent it is important to note that the tensor Djkin can 
depend on external force fields. This implies that the 
eigenvalue uij of the matrix W can depend on the fric- 
tion coefficient a. If we were to assume the quantity 
ujj to be independent of the friction coefficient a, then 
we obtain the expression of the Lyapunov exponents as 

\f ] = ±\u J /{2m)\ 2 / 3 -a/2 + 0{a 2 ) from Eq. (||) in the 

case of \a\ < -\/3 \ujj/m\ 2 ^ 3 . However this is not consis- 
tent with the numerical results in a deterministic many- 
hard-disk system with a color field in which the negative 
Lyapunov exponents rather increase as the value of the 
friction coefficient increases Jl2) . This consideration sug- 
gests that we should not neglect the external force field 
dependence of the correlation amplitude DjUn at least in 
the color field case. The dependence of the tensor Djkin 
on the shear rate and the external force fields should be 
a subject for a separated paper, although the conjugate 
pairing rule of the Lyapunov spectrum is correct regard- 
less of their dependence as shown in this paper. 



VI. CONCLUSION AND REMARKS 

In this paper we have applied the master equation 
approach to Lyapunov spectra to non-equilibrium iso- 
kinetic thermostatted systems in order to discuss a conju- 
gate pairing rule. We considered two-dimensional many- 
particle system with Gaussian white random interactions 
between the particles. In this system the positive Lya- 
punov exponents are calculated by a (forward) Fokker- 
Planck equation for the tangent vector dynamics. Wc 
proposed a method to calculate the negative Lyapunov 
exponents by a time-reversed master equation, especially 
the anti-Fokker-Planck equation where the diffusion term 
has the opposite sign to the forward Fokker-Planck equa- 
tion. Using the Lyapunov exponents calculated by these 
two Fokker-Planck equations we showed the conjugate 
pairing rule of the Lyapunov spectrum for iso-kinetic 
thermostatted systems with a shear field given by the 
Sllod equation in the thermodynamic limit. We also 
gave an concrete form to connect the Lyapunov expo- 
nents with the time-correlation of the interaction matrix 
in a thermostatted system without a shear field. 

We discussed the conjugate pairing rule based on the 
iso-kinetic thermostat in the thermodynamic limit. How- 
ever it is known that the iso-kinetic thermostat is for- 
mally equivalent to other thermostats such as the iso- 
energetic thermostat in the thermodynamic limit [p5| . 
In this sense our result should be correct in systems with 
such other thermostats, more explicitly as far as the fric- 
tion coefficient can be regarded as a constant even in a 
finite number of particle systems. 

In order to get the ani- Fokker-Planck equation we used 
the fact that the shear rate 7 changes its sign in the time 
reversed motion. However this time-reversed change of 
the sign of the shear rate to get the anti-Fokker-Planck 
equation may not be essential to obtain the negative Lya- 
punov exponents, if the Lyapunov exponents are even 
functions of the shear rate. We have not proved the in- 



variance of the Lyapunov exponents under the transfor- 
mation 7 — > —7 in this paper, but the invariance of Eq. 
( p2| ) under the transformation Q — ► — O implies that the 
Lyapunov exponents are invariant under this transforma- 
tion. 

We can show that all the Lyapunov exponents Aj + ' 

(A^ ') are non-negative (non-positive) in the case of 
7 = (See Appendix ^.). This implies that in this case 
the number of the positive Lyapunov exponents should 
be equal to the number of the negative Lyapunov expo- 
nents. However we have not proved that it is also correct 
in the presence of a shear field: 7 7^ 0. Concerning this 
point we should notice that a numerical calculation of the 
Lyapunov spectrum for the Sllod equation (H) and (|^) 
showed that in the case of a high shear rate the number 
of the positive Lyapunov exponents can be less than the 
number of the negative Lyapunov exponents |30| . There- 
fore it should be interesting to check whether the master 
equation approach to the Lyapunov spectrum can de- 
scribe such a situation or not. 

It should be noted that the discussion of this paper 
does not conclude that the conjugate pairing rule of the 
Lyapunov spectrum for the thermostatted system with a 
shear field must be satisfied rigorously in deterministic 
chaotic systems. To show the conjugate pairing rule in 
this paper we assumed the Gaussian white randomness 
( |l0| ) and ( [Tl| ) for the particle interactions, and there is 
no guarantee that we can justify the conjugate pairing 
rule by the master equation approach under a more gen- 
eral random interaction of particles, especially under the 
randomness with a long time correlation which leads to 
a more general master equation for the tangent vector 
rather than a simple Fokker-Planck equation. A gen- 
eralization of the conjugate pairing rule by the master 
equation approach to a more general random interaction 
remains as one of the important future problems. 
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APPENDIX A: MASTER EQUATION FOR THE 
TANGENT VECTOR 

In this appendix we derive the Fokker-Planck equation 
( |l2| ) for the tangent vector space. Using the Kramers- 
Moyal expansion the dynamics of the probability density 
p(+>(Sr,t) is given by §7) 

00 2N 2N 2N 



dt 



where • (£T,t) is defined by 



= £££•••£(-1)" 

n=l 31=1,73=1 j n =l 



d5r h dsr n ■■■dsr jn 



(A.l) 



6 



(6T,t) 



^\im\{[5r h {t + s)-5r h {t)\ 



nl s^o s 



x[sr h (t + s)-sr h (t)} 
■■■x[sr jri (t + s)-6r jn (t)}) 



(A.2) 



and STj(t) is the j-th component of the tangent vector 
ST(t). 

Using Eq. (|5|) we obtain 



ST(t 



6T(t) 





■ i-t+s 


|r exp 


J dr C(t) 



E 

n=l 



dr n J dr n -i ■ ■ 



1 ST(t) 



dri 



xC(T n )£(T n - 1 )'--C(T 1 )5T(t). (A.3) 
It follows from Eqs. (@), ©, @, jA~2|) and ( |A~3| ) that 



lim- ([JlYt + s) -ST(t)]) 



5r(t)=<5r 



lim — 



dr (£(r)> <Sr 



/ <Mq + Sp/m 
\ tySp/m 



(A.4) 



and 



sP)(ir,t) = (sJ'(ir,()) 

= limi/[5r(t + s )-<5r(t)] 

x [<sr(t + «) - ar(t)]' J 

1 z"^ 5 /■*+* 
= lim — / (i/t / 
2s J t J t 

x //:( K ) ( 5r ( 5r T /:(T) :! 



sr(t)=sr 



02JV 02AT 

2 at r?(5q) 



where f?(<5q) = (77^ (dq)) is defined by 



N N 



(A.5) 



(A.6) 



1 V—\ 



Here the only non-ze ro c ontributions come from the 
n = 1 term of Eq. (|A.3| ). For general n, the num- 
ber of delta functions from Eq. (^Tj) must be only one 
less than the number of time integrals, to give a non- 
zero contribution. It is straightforward to show that this 



never happens for n > 1. Concerning the terms including 
Bj"l»—.7 (ST, i), n = 3, 4, • • • in the right-hand side of Eq. 
( A.l) we obtain 



S S 2 -i n ( <Jr '*)=° ! n = 3,4,-. 



(A.7) 



because of the Gaussian white properti es (Eo[) and ( 0) o f 
the r ando m matrix R(t). Using Eqs. ( |A.l| ), (A. 4), ( |A .5|) 
and (A.7) we obtain the Fokker-Planck equation Q12|). 



APPENDIX B: EQUATION FOR THE MATRIX 

f (±) 

In this appendix we give details of the derivation of 
Eq. @ from Eqs. @, @ and (jl|). We start this 
derivation by introducing the N x N matrices T^it) 
and Q^(t) defined by 



T^(t)^(S q Sp T ) 



T\(±) 
t 



g(±)(t)^(5p6 P T )f. 

Equations (|l2|) and rtlS) lead to the equations 



(B.l) 
(B.2) 



1 r 



(B.3) 



dT^\t) 



dt 



m 



(B.4) 



*ac±) (<) + g(±) ± p |t(±) (t) } 



(B.5) 



for the matrices T^i), J" (±) (t) and £ (±) (i) with the 
oper ator 2? defined by Eq. (p4|). Here, to derive Eq. 
(B.5) we used the relation 



D 



Injk -Djkln 



(B.6) 



which is derived from the definition (|Tl|) of the tensor 
Djki n a nd th e sym metry prop erty of the matrix R(t). 
Eqs. fO|), (IO) and (R5) are equivalent to 



df(^(t) 



dt 



= T (±) [t)P{t) T + P{t)T {±) (t) T , (B.7) 



dt 



(B. 
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dt 



= ±V t 



{f(±)(t)} 



(B.9) 



where f (±)(t), and g (±) (*) are defined by 

fWfi) = e-*!^)^-*'', (B.10) 



^ ± )(t) = e -* t ^ ± )(f)e-* T< 



^(i) = e-* f 5( ± )(*)e-* r, 1 
and P(t) is the (2iV) x (27V) matrix defined by 



(B.ff) 
(B.12) 

(B.13) 



with the matrix O defined by Eq. (p3|), and D t {- ■ ■} is 
defined by 

V t {X} = e- m V {e* ! Ie* Tt } e -* T * (B.14) 

for any (2N) x (2N) matrix X. Here we used the relation 

= (B.15) 

so that we have the equation exp{— <I>t} ■ exp{^t} = 
exp{-($ - *)*}. 

Noting that the matrix Q^\t) is symmetric and the 
inverse matrix of the matrix P(t) is given by P{t)~ l = 
m,exp{f2t}, we obtain 



2Q^{t) = P(ty^-^l + ^ [P(t)-n ' 

^d 2 T(±)(t) r , . , 



p(ty 



dt 2 



dt 

[pitr 1 ]- 



+P(t)- 1 f<- ± \t)Q T 

+o^±)(t) T [Pity 1 ] 7 * 



(B.16) 



by us ing E qs 



and (B.13). Besides, using 
Eqs. ( jESTD , (|B.8|) and ( gig ) we obtain 

^{p(t)- 1 ^ ± )(t)n T + ^( ± )(t) T [p(f)- x ] T } 



+g( ± >(t)f> T + ^ ± )(t) 



(B.17) 



where we again used the relation (/ (jjj (t ) T = Q^(t). It 
follows from Eqs. (|R7|), |tf ) and (^J) that 



±22?. 



{t (±) w} 



dt 



-QP(t) 



dt 2 
_ 1 df( ± )(t) 



[PW^O 3 



±5 (±) (t)r! T + ^ (±) (t). 



(B.18) 



Takin g the time-differ entia l of both the sides of Eq. 
( B.18 ), and using Eq. ( |B.9| ) we obtain 



J2 J2f(±)( f \ 



dt 2 

d D ^-id* (±) (*) 



- Q dt P ^ 



fit 



[P(t)- 1 ] T ^ 



±2^{f( ± )(t ) } 

±2? t {f (±)(t)} T ± Q2? t {f (±)(t)| 



(B.19) 



This is the equation for the quantity T^^t) only. 

Now we derive the equation for T^(t) defined by Eq 
(0) from Eq. ( |B.1S| ) for fW(t) defined by Eq. QB.10[) 
We note 



= 2$ + a/ 2 Ar 
= -2* - aI 2 N 



(B.20) 
(B.21) 



which is derived from Eqs. (@), (|) and (||. Using 
Eqs. (||), flB.lCj ) and JB.20D the matrices T^(t) are 
connected with the matrix T^(t) by 



fW(t) = e- ot / 2 f( ± )(±t) e -^ 



t/2 



(B.22) 



We introduce the multiplication X (g)Y oi X and Y which 
is defined by 



X $ Y = 1 [jf Y + (XY) T ] 



(B.23) 



for any square matrices X and Y of the same size. This 
multiplication is used in the relation 



dt y ' 



,±nt/2 



dZ(t) 
dt 



± O <g> Z(t) 



,±n T t/2 



(B.24) 



satisfied by any (2iV) x (2N ) sym metri c mat rix Z( t) as a 
function oft. Noting Eqs. ( |B.23j ) and ( |B.24| ), Eq. ( |B.22| ) 
leads to 



P(t) 



_ x fif (±)(t) 



fit 



[Pity 1 ] 



m 2 e nt ' 2 



^M-^f^(±t) 
at 



e° T */ 2 , (B.25) 



P(t) 



</t 2 



P(t)" 



f/ 2 f( ± )(±t) fif( ± )(±t) 



dt 2 
+0c 



fit 



n® f (±) (±t)]}e" T * /2 (B.26) 
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where we used the relation T^ ± \t) T = T^(t). More- 
over T> t operated on the matrix fW(i) is connected with 
T> operated on the matrix T^(±t) as 

V t {t (±) (t)} = e nt / 2 V {?&(±tj\ e " 4 / 2 (B.27) 



where we use d Eqs . (PU|I^) fjB.14|) and ( |B.21 
serting Eqs. (|B 25|) pl26j ) and ( ^27| ) into Eq. 
and using Eq. ( |B.24| ) we obtain 



In- 



B.1E) 



d 4 f(±)(±0 



+0 



2si< 



fit 



d 2 f ( - ± \±t) 



dt 2 



d 2 f^(±t) nT 



dt 2 



9. 



SJ®f (±) (it) 



si 



si<g> f (±) (±t) 

' dt (±) (±t) | 



si T 
o. 



(B.28) 



Eq. (B.28) is equivalent to 
d 4 f (±)(±i) 



dt 4 
1 

~2 
1 

f 16 



= 0. 



2 d 2 t ( -^(±t) t d 2 f^(±t) 



dt 2 



dt 2 



si 2 



S1 4 T (±) (±i) - 2S1 2 T (±) (±i) 



-T (±) (±i) (SI 4 



dt 



(B.29) 



By exchanging i with ±t in Eq. (B.29) we obtain Eq 



APPENDIX C: LYAPUNOV EXPONENTS IN 
THE COLOR FIELD CASE 

In this appendix we consider the case with no shear 
field using the condition ( p9| ) , and derive Eq. ( |3l| ) under 
the assumption ([50|). We also give a derivation of Eq. 
(H) briefly. 

Under the condition ( |29| ) , the matrix SI defined by Eq. 
( p3| ) is simply an identical matrix multiplied by a con- 
stant and is given by 

SI = aI 2Nl (C.l) 

and the matrix P(t) defined by Eq. ( B.13| ) and the op- 



erator T> t {- • •} defined by Eq. (B.14) are given by 
P(t) = -e- &t I 2N , 



(C.2) 



i> t {X} = e 2Blt V{X} (C.3) 
for any (2N) x (2iV) matrix X. Noting Eqs. dcj) and 



(C.3 
(B 



and the relation f^(t) = T^(t), Eqs. (| 
and flB.9|) are simply attributed into 



dT(±)(t) 



1 



dt 



df^{t) 



at m 



(C.4) 



(C.5) 



dQ {±) {t) , 
dt 



Noti ng tha t the matrix Q^\t) is symmetric, Eqs. (C.4), 
(|C.5|) and (|C.6|) lead to the differential equation 



da t d at d r(±) = ± 2 e2at p ( T (±) W \ (c.7) 
ai eu eft rrr I > 

for the matrix T^'^t) only. 

Now we consider the derivation of the equation for the 
matrix f^(t) (defined by Eq. ©) from Eq. (|07[) fo r 
the matrix T^t). It follows from Eqs. ((2lJ) and ( C?fy 
that 



(C.8) 



By exchanging t with ±i in Eq. ( p.8[ ) and using Eqs. 
@ and © we obtain Eq. @. 

Using Eqs. @ and (§f) the real function ^ (±) (t) of f 
defined by 



,±at 



satisfies the equation 

,(±) 



(C.9) 



d^f>(t) _ 3 d£f>(t) 2w? (±) 



A 3 



(ft 



_i^(t) = 0. (CIO) 



The real solution of the linear differential equation ( C.10| ) 
is expressed as 



^\t) = J2Re{v^e^} (C.ll) 



fe=i 



where j — 1,2,3 are constants determined by the 

initial condition, and Q , j — 1>2,3 are the three solu- 
tions of the equation 

„ , 2w 2 
C 3 - a 2 C - -f = 



(C.12) 
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for £. Here Re{AT} means to take the real part of 
any imaginary number X. We sort the quantities Q, 

j = 1,2,3 as Reld 15 } > Re{Cf } } > Re{Cf } }, so that 
using Eqs. (|27|), ( |3^ ) and (C.S) the Lyapunov exponent 
A^ is expressed as 



Af > = ± 



lim — In — 



(±: 



(t)e 



=pat 



d ±} (o) 



(C.13) 



It follows from Eq. ( gjg ) that C = <f } 

is a real solu- 



tion of Eq. ( C.12j ) and satishes the conditions (. 



(i) 



(i) 



A,- 



a 

3A~ 



(C.14) 



We can 



with the quantity Aj defined by Eq. 
check that in the case of \a\ < \ojj/m\ 2 / 3 the quan- 
tities Aj and ^j 1 ^ are both real numbers, and in the 
case of I a I > v3 \ujj/m\ 2 / 3 the quantity Aj can be an 
imaginary number but the quantity given by Eq. 



(C.14) is still a real number and satisfies the condition 



linx. 



A 1 



= |a|. Using Eqs. ( |C.13|) and <\C.l4) wc 



obtain Eq. (64). 
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